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Abstract: The rupture of thin-cap fibroatheroma accounts for most acute 
coronary events. Optical Coherence Tomography (OCT) allows 
quantification of fibrous cap (FC) thickness in vivo. Conventional manual 
analysis, by visually determining the thinnest part of the FC is subject to 
inter-observer variability and does not capture the 3-D morphology of the 
FC. We propose and validate a computer-aided method that allows 
volumetric analysis of FC. The radial FC boundary is semi-automatically 
segmented using a dynamic programming algorithm. The thickness at every 
point of the FC boundary, along with 3-D morphology of the FC, can be 
quantified. The method was validated against three experienced OCT image 
analysts in 14 lipid-rich lesions. The proposed method may advance our 
understanding of the mechanisms behind plaque rupture and improve 
disease management. 
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1. Introduction 

Coronary artery disease (CAD) is the most frequent cause of death worldwide. Each year, an 
estimate of 785,000 Americans will have a new coronary attack, and about 470,000 will have 
a recurrent event [1]. Most acute coronary events result from rupture of the protective fibrous 
cap (FC) overlying an atherosclerotic plaque [2]. Thin cap fibroatheroma (TCFA) has been 
identified as the most frequent precursor lesion that leads to plaque rupture [3]. TCFA has 
been characterized as a plaque containing a large lipid necrotic core covered by a thin FC 
infiltrated by macrophages [3,4]. An ex vivo morphometric assessment of 41 ruptured 
coronary plaques revealed that the mean cap thickness was 23 + 19 |Lim and 95% of these FCs 
measured less than 64 |im [5]. Because of this, a thickness of 65|im has been considered the 
instability threshold, and thickness of the FC has been considered one of the major 
morphometric determinants of those plaques prone to rupture [3,4]. Identifying these so-called 
"vulnerable plaques" before symptoms arise, may allow the adoption of preventive 
therapeutic measures to avoid subsequent myocardial infarctions and sudden deaths. This 
largely relies on improvement of imaging technologies that allow better characterization of 
fibroatheroma in vivo. 

Intravascular Optical Coherence Tomography (OCT) is a high resolution (10-20 jim) 
imaging modality that can image the depth-resolved profile of coronary artery to about 1- 
2mm [6,7]. OCT has better sensitivity and specificity for detection of lipid plaques as 
compared to intravascular ultrasound (IVUS) [8] and is currently the only imaging modality 
that can measure the FC thickness in vivo [9]. OCT has already been used in several in-vivo 
studies to assess the correlation between the incidence of TCFA/FC thickness and clinical 
presentations in living patients [10-17]. 

The current accepted standard method for assessing FC thickness using OCT images is 
based on single measurements of the thinnest portion of the FC [9,18]. In practice, the 
segment of interest of the coronary vessel is imaged in 3D using a catheter probe and a pull- 
back procedure. The image data are subsequently analyzed manually by expert analysts. The 
cross-sectional image within a field of view encompassing a lipid-rich plaque where the FC 
thickness appears to be thinnest is visually identified. The OCT criteria for identification of a 
lipid-rich plaque are the presence of a signal-poor region with diffuse boundaries and high 
attenuation [19]. In the selected cross section, the thinnest portion of the FC is visually 
determined and its thickness is manually measured. In order to decrease variability, averaging 
multiple FC thickness measurements has been proposed [13,20]. 

However, there are two limitations associated with this method. First, the measured FC 
thickness is subject to inter-observer variability due to visual uncertainties involved in the 
entire process. Second, the FC is intrinsically a 3-D structure and single pointZ-single frame 
measurements are unable to characterize the volumetric nature of FC. Although the thickness 
of FC has been suggested a risk factor for plaque rupture, the exact mechanisms behind FC 
thinning and plaque rupture are unknown. A more consistent and comprehensive method for 
FC quantification may better elucidate the unknown factors associated with plaque rupture. In 
this paper, we present a new, computer-aided algorithm for semi-automatic volumetric 
quantification of the FC. With this method, the circumferential distribution of the lipid 
plaques/FC was first identified by human analysts in all frames containing FC, then the FC 
boundaries in the radial dimensions were automatically segmented and volumetric 
quantification was performed on the segmented FC. In the following sections, we first 
describe how the FC boundaries can be accurately and robustly segmented using a dynamic 
progranmiing algorithm. We validate the method by comparing its accuracy and variability to 
that provided by three experienced OCT analysts. To demonstrate the significance of the 
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method, we assess the intra-observer and inter-observer variabiHty of the conventional single 
measurement by human analysts. Lastly, we introduce quantitative volumetric metrics that 
can be derived from the segmentation results. 

2. Methods 

2.1. Population, image acquisition and study protocol 

The population of this study is comprised of 10 elective patients who underwent OCT 
evaluation of target coronary arteries prior to percutaneous coronary intervention, between 
September and October 2010. All OCT images were acquired by a commercial Fourier 
Domain OCT system (C7-XRTM OCT Intravascular Imaging System, St. Jude Medical, St. 
Paul, Minnesota) at University Hospitals at Case Medical Center (Cleveland, OH), digitally 
stored and de-identified by the Institution's Cardiovascular Imaging Core Laboratory. The 
axial resolution of the system is about 15 jim. The scan characteristics of the system are: 
50000 lines/s, 504 lines/frame, yielding 100 fps and 20 nmi/s puUback speed yielding a 200 
|Lim frame interval. Lipid-rich plaques were selected by an interventional cardiologist 
experienced in reading OCT images according to the presence of signal-poor regions with 
diffuse borders under OCT [19]. Calcified plaques may also produce a bright interface 
overlying a signal-poor region which, opposed to lipid-rich plaques, has a heterogeneous 
pattern with well-defined sharp boundaries [19], and were excluded from this study. Once the 
lipid-rich plaque was identified, all of the consecutive frames in the puUback containing the 
lesion were selected. A total of 323 images from 14 lipid rich lesions were included and used 
for validation. 

2.2. Image analysis method 

The OCT raw data were logarithmically compressed and transformed from the polar 
coordinates (9, r) to the Cartesian coordinates (x, y) consisting of 1024 by 1024 pixels (pixel 
size: 9.4 by 9.4 |im). The major steps for volumetric quantification of FC are (1) segmentation 
of FC boundaries in all frames containing the lesion, and (2) quantification based on the 
segmentation. 

FC is delineated by a luminal boundary and an abluminal boundary. In the non-ruptured 
plaque the FC protects the underlying lipid/necrotic core from contacting the circulating 
blood. Therefore, the luminal boundary coincides with the vessel lumen contour. The 
abluminal boundary is conmionly described as the "diffuse border" created by the interface 
between the FC and the underlying lipid pool. We adopt a sequential approach to first detect 
the vessel luminal boundary, based on which we can then track the FC abluminal boundary. 
The lumen contour will also be used for validation and quantification purposes. 

Searching for the vessel luminal boundary and FC abluminal boundary in the polar 
coordinates of OCT images fit the graph search family of optimization problems naturally. 
The lumen boundary is unique and the accumulated optical intensity difference between the 
pixels from the vessel and luminal side along the contour is maximum. The FC abluminal 
boundary is the optimum boundary that best separates the FC from the underlying lipid 
plaque. 

2.2.1. Dynamic programming 

Dynamic programming (DP) is a general technique used to solve optimization problems [21]. 
The basic concept is to find the global optimal solution to the original problem by building on 
optimum solutions to subproblems. It is very robust in the presence of noise and this is 
attractive to OCT image analysis because OCT images often suffer from speckle noise and 
various artifacts. DP has been used previously in intravascular OCT to correct non-uniform 
rotation distortion [22]. Consider the OCT images in polar coordinates (9, r) where 9 is angle 
and r is depth. We assign each pixel at row / and colunm j in the OCT image an objective 
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function favoring the characteristics of lumen or FC abluminal boundary. Given the 
starting row (angle) and ending row 6^^ of the boundary, we are essentially searching a 
path Pq^^Oj^ from 6^ to 6^^ with the optimal cost C. We can break the problem into 
subproblems such that Pq^^Oj^ is always coming from Pq^_^q^_^ with some connectivity 
constraint. Therefore, we have the following recursive function: 

c{ij)= max [c{i-i,f)+f{ij)] e^<i<e^ 

i-n<i*<j+n 

where C(/, j) is the accumulated cost from row 6^ to point (/, 7), 7* is adjacent to j, and n 
specifies connectivity. The global optimum boundary can be found by selecting the point in 
row where the accumulated cost is maximum and back tracking the path. Equation (1) 
defines the common method used for lumen and FC abluminal boundary segmentation. Only 
the objective function is different for different tasks. 

2.2.2. Lumen and guide wire segmentation 

The objective function for lumen segmentation can be defined straightforwardly as the pixel 
value difference between the vessel side and the luminal side: 

f_lumen(i, j) = T(i, j < j + m^)- T(i, j-m^K < j) (2) 

where / refers to the average of pixel value, and and are the length of windows on the 

vessel and luminal side, respectively. In this study, we chose -m^ =0.1 nmi. Under this 

definition, pixels of the guide wire reflections may obscure the lumen boundary and need to 
be excluded from lumen segmentation. We have proposed a guide wire segmentation method 
that is able to extract the global optimum guide wire positions of all the frames in the entire 
puUback at once [23]. Briefly, a 2D longitudinal image of average A line values plotted as a 
function of frame number and 6 is created (Fig. 1(b)). The regions of guide wire shadow 
become a continuous dark band. Similar to lumen segmentation, an objective function of pixel 
value difference is applied to the two boundaries of the dark band, but with different signs. 
DP is then applied twice to find the two boundaries. 

2.2.3. FC abluminal boundary determination 

Although OCT has been suggested to have the ability to quantify FC thickness accurately 
[18], a quantitative definition of the FC abluminal boundary under OCT has not been 
proposed. The FC abluminal boundary is usually characterized by a gradual transition of pixel 
intensity from bright to dark, and visual determination of the optimal transition point can be 
cumbersome, inaccurate and subject to high variability. On the other hand, given an objective 
function favoring the appropriate properties of FC, computers can objectively identify the 
location of the abluminal boundary more consistently than humans. The FC has been defined 
histologically as a distinct layer of connective tissue covering the lipid core, and consists 
purely of smooth muscle cells in a collagenous-proteoglycan matrix, with varying degrees of 
infiltration by macrophages and lymphocytes [4]. The fibrous tissue appears bright by OCT 
and has a low attenuation coefficient, whereas the lipid pool appears dark and strongly 
attenuates the light [24,25]. Therefore, we suggest that the FC abluminal boundary has a high 
intensity difference between the FC and lipid pool, and also a high gradient. Hence, we define 
the following objective function: 

f_FC{U j) = T{U j-d,< < j) - T{U j <j.<d^)- Aju (3) 
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Fig. 1. Guide wire and lumen segmentation, (a) OCT image in polar coordinates. The guide 
wire has a high reflection followed by a long dark shadow. The average intensity of each A 
hne was computed and the 2D image became a single Hne illustrated as the dotted line in (b). 
(b) The guide wire region became a continuous dark band traversing the whole puUback after 
the average intensity of each frame was combined into a single longitudinal image. Dynamic 
programming was performed on this image to extract the two boundaries of the guide wire 
shadow, (c) Corresponding anatomical view OCT image of (a), (d) Lumen segmentation using 
the dynamic programming algorithm. 

where and d^^^ are predefined depths to calculate pixel intensity difference, ju is the slope 

of pixel value attenuation extracted from the A line segment of length L across (i,j), and /I is a 
weighting term. The parameter values d^ =15 |am, d^^ =0.38 mm, /I = 7, and L = 38 jim 

used in the validation set were determined experimentally using a separate training data set 
(n = 3) different from the 14 lesions selected for validation. The ROI for analysis is bounded 
in the radial dimension by the previously segmented lumen boundary and a predefined 
maximum depth d^^^ , and in the circumferential dimension by a user-defined starting and 
ending angle encompassing the selected lesion. 

2.2.4. Volumetric quantification 

With the fully segmented FC, we can quantify the thickness at each point of the FC luminal 
boundary, defined as the minimum distance from this point to the FC abluminal boundary. 
The conventional metric, the minimum FC thickness (MCT), of a single lesion can be simply 
found by searching for the minimum thickness out of all the points on the FC luminal 
boundary in all consecutive frames covering this lesion. The mean FC thickness (MeanCT) of 
a lesion is the average thickness of all points on the FC luminal boundary. 

In addition, volumetric metrics of the FC can be defined as follows. The FC surface area 
(SA) of a lesion can be calculated as the product of the frame interval and the arc length of FC 
summed over involved frames. The arc length can be determined from the radius of FC 
luminal boundary with reference to the centroid of the lumen. The absolute dind fractional FC 
categorical surface area (ACSA and FCSA) of a lesion can be calculated as the absolute and 
relative FC area in a thickness category. In this study, we classified FC thickness into 3 
categories: <65 jim, 65-150 jim and >150 jim. This definition is based on both pathology 
studies and in vivo clinical studies using OCT [4,17]. However, further evaluation of the 
significance of these cut-off values will be important in the future. FC volume (Vol) of a 
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lesion is calculated from cross sectional FC areas in individual frames using Simpson's rule. 
FC volume density (VD) of a lesion is the FC volume normalized by the lesion length. FC 
fractional luminal area (FLA) is defined as the percentage of luminal area occupied by the 
FC. Categorical FC fractional luminal area (CFLA) is defined in a similar way but only 
counting the FC area in a thickness category. Cross correlation between the volumetric 
metrics will be used to evaluate redundancy (collinearity). 

2.3. Validation experiment 

Three experienced OCT analysts were involved in the validation experiments. A graphical 
user interface written in MATLAB (MathWorks, Inc) and C++ for visualizing and analyzing 
OCT images was developed and provided to the analysts. Importantly, proper calibration of 
all OCT image runs through adjustment of the Z-offset [6] was performed by a single 
operator, and kept the same by the other two analysts and by the computer algorithm. 
Likewise, the same operator was responsible for selecting all lipid-rich lesions included in the 
analysis. 

In order to allow a comprehensive understanding of the importance of a full volumetric 
segmentation of the FC, inter-observer variability of FC thickness quantification of all 14 
lesions were initially determined by the single-frame, single -measurement concept among the 
three analysts. In this process, each analyst independently selected the image frame where 
they suspected the FC was thinnest and determined the FC thickness by a single measurement. 

Subsequently, all 323 cross-sectional images from the same 14 lipid-rich lesions were 
manually segmented by analyst 1 and by the computer algorithm. The accuracy of the 
computer algorithm was evaluated by comparison with the manual segmentation. In order to 
assess intra-observer variability, all images were re-analyzed by operator 1 two weeks later. 
To assess inter-observer variability of manual segmentation of FC boundaries, 50 randomly 
selected cross-sections were analyzed by operators 2 and 3. At all times, each analyst was 
blinded to the analysis results performed by the other analysts and the computer algorithm. 
Only the common region of FC selected by all three analysts was used for comparison. If the 
FC was blocked by the guide wire shadow, only the shadowed region was excluded for 
validation. 

Since the algorithm determines both the luminal and abluminal boundaries of FC, we 
evaluated the accuracy of every point on the two boundaries as illustrated in Fig. 2. If we 
project n rays from the lumen centroid to the FC at one degree intervals, and define the 
distance from the lumen centroid to the FC luminal and abluminal boundary along ray / 
as d_LB. and d_ALB. , respectively, we can use the mean absolute difference (MAD) between 

these distances determined from the semi-automatic method and manual segmentation to 
assess the agreement: 




Fig. 2. Methodology for validation of fibrous cap luminal and abluminal boundary 
determination. The distance from the two boundaries to the lumen centroid was defined as 

d_LB. and d_ALB. , respectively, and compared between the computer algorithm and human 

analysts. This distance was calculated along rays at intervals of one degree. 
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n i=\ 



d_LB-d_LB. 

Manual Auto 



MAD,,,=-T 

n i=\ 



d_ALB-d_ALB. 

Manual Auto 



(4) 



The mean signed differences (MSD) can be similarly defined. This stringent metric allows 
complete comparison of every point on the FC boundary. 

3. Results 

3.1. Intra- and inter-observer variability of human analysts 
3.1.1. Single-point measurement of minimum cap thickness 

The intra-ob server variability in determining the minimum cap thickness (MCT) was 14.5 ± 
11.7 jim for analyst 1. Figure 3 illustrates the manually measured MCT by the three analysts. 
The concordance correlation coefficient (CCC) of the frame selection was excellent (ccc = 
0.9958) among the three analysts, with a mean absolute difference of the selected frame for 
assessing the MCT of 3.1 ±4.1 frames. However, the mean absolute difference between the 
three analysts in assessing the MCT of the 14 lesions was 12.8 ± 8.2 jim, with a very poor 
CCC (ccc = 0.2450). Of note, when the cut-off point of 65 jam was used to define a TCFA, 
agreement among all three analysts was reached in only four lesions and at least one analyst 
held a different assessment in 10 out of the 14 lesions. 



CCC=0.9958 



1 . « t X 



■ Analyst 1 
^ Analyst 2 
• Analysts 



Lesion Number 

(a) 



Q. 65 
CO 
O 



CCC=0.2450 



Lesion Number 
(b) 



Fig. 3. Manual measurements of the minimum cap thickness of the 14 lesions by the three 
analysts, (a): The difference of the selected frame for measuring the minimum cap thickness by 
the three analysts was plotted vs. the lesion number. The frame difference is against analyst 1. 
Frame difference of 0 indicates no difference, (b): The minimum cap thickness measured on 
the selected frame in (a). CCC: concordance correlation coefficient. CCC<0.4 is considered to 
be poor agreement. Dashed line: conventional cut-off value of 65 |am for definition of TCFA. 
Under this definition, one of the three analysts held a different assessment of TCFA for 10 out 
of 14 lesions. 

3.1.2. Manual segmentation of FC 

The Mean Absolute Difference (MAD) of the segmentation performed by human analyst 1 in 
two different time points was 20.7 ± 11.3 jam. The MAD between the three human analysts 
was 30.3 ± 27.3 |im. The angular difference of the selected circumferential distribution of FC 
for the three analysts was 9.6 ± 9.9 degrees for all lesions except for one outlier (63.0 
degrees). 

3.2. Performance of the algorithm 

The average processing time for lumen segmentation, guide wire segmentation and FC 
abluminal boundary determination in a single frame was 0.4s, 0.07s and 0.2-0.4s, respectively 
(with a duo core 2.00GHz CPU). 
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The MAD of the FC thickness between the algorithm and analyst 1 in assessing all the 323 
images was 27.3 ± 26.7 jam. The MAD of the luminal and abluminal boundary was 15.7 ± 
23.4 |im and 25.3 ± 31.4 |im, respectively. The Mean Signed Difference (MSD) of the FC 
thickness, luminal and abluminal boundary between the algorithm and analyst 1 was 2.9 ± 
38.0 |am, 6.0 ± 27.5 [xm and -3.2 ± 40.2 |im, respectively. 

Figure 4 shows the MAD and MSD comparison of the localization of the luminal and 
abluminal boundaries and the thickness of the FC by the algorithm and by the three analysts. 
The MAD between the algorithm and analysts was comparable to that between analysts. The 
MAD for luminal boundary was less than that for the abluminal boundary. From the MSD 
results in Fig. 4(b), the bias is around one pixel or less, except between analyst 3 and the 
algorithm and between analyst 3 and the other two analysts in determining the FC abluminal 
boundary (and consequently the FC thickness). 
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Fig. 4. (a): Mean absolute difference (MAD) between automated and manual measurements by 
three observers, (b): Mean signed difference (MSD). Left: fibrous cap thickness determination. 
Middle: Luminal boundary determination. Right: Abluminal boundary determination. O: 
observer. Auto: computer algorithm. 

Figure 5 illustrates a typical example where the semi-automatically determined contour 
was very similar to all three analysts. While the analysts tended to smooth the FC boundary, 
the algorithm exactly followed the border where the intensity difference and gradient were 
high. Figure 6 shows an example where the algorithm shows disagreement with the three 
analysts. Compared to the algorithm determined segmentation, the FC traced by analyst 2 is 
significantly thinner, whereas the one traced by analyst 3 around 9-10 o'clock is thicker. 
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Fig. 5. In this case, the semi-automatically segmented fibrous cap boundary showed good 
agreement with all three human analysts. 




Fig. 6. In this case, the algorithm shows disagreement with the 3 analysts. The traced FC from 
analyst 2 is significantly thinner than that from the algorithm. In comparison, the FC traced by 
analyst 3 around 9-10 o'clock is significantly thicker. Analyst 1 selected more circumferential 
distribution of FC for segmentation. 

3.3. Volumetric metrics vs. minimum cap thickness 

Among the 14 lesions, there was low correlation between all the volumetric metrics and the 
conventional metric MCT (Table 1). This indicates that these volumetric metrics characterize 
independent morphological information about FC. Between the volumetric metrics, there was 
moderate correlation between ACS A (<65 |im) and MeanCT (R = -0.67, p<0.01), between 
FCSA (<65 |im) and MeanCT (R = -0.76, p<0.01), between CFLA (<65 |im) and MeanCT 
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(R = -0.65, p<0.01), between ACSA (<65 |im) and SA(R = 0.74, p<0.01), between ACSA 
(<65 |Lim) and Vol (R = 0.64, p<0.01), between FLA (<65 |im) and SA (R = 0.64, p<0.01), 
between VD and SA (R = 0.77, p<0.01), between VD and FLA (R = 0.77, p<0.001), and 
between VD and Vol (R = 0.82, p<0.001). There was strong correlation between Vol and SA 
(r = 0.99,p<0.0001). 

3.4. Three dimensional visualization ofFC 

The segmented FC can be visualized in 3D with a continuous colormap, from blue to green to 
red, indicating the FC thickness ranging from 300 |Lim to one pixel length. Figure 7 illustrates 
2 cases, both with TCFA, with FC thickness rendered in 3D. If assessed using only the 
conventional methodology, the morphological differences between the cases would not be 
apparent. However, the 3-D visualization demonstrates dramatically different characteristics. 
In particular, the vessel shown in the lower panel has a significantly higher fraction of thin FC 
as compared to the one in the upper panel. 

Table 1. Correlation coefficients between the volumetric metrics in 14 lesions'" 





MCT 


MeanCT 


ACSA 
(<65 |Lim) 


FCSA 
(<65 |Lim) 


SA 


FLA 


CFLA 
(<65 |Lim) 


Vol 


VD 


0.02 


0.18 


0.35 


-0.49 


0.77* 


0.77** 


-0.26 


0.82** 


Vol 


-0.31 


-0.14 


0.64* 


-0.30 




0.63 


0.08 




CFLA 
(<65 jLim) 


-0.23 


-0.65* 


0.39 


0.47 


0.17 


0.11 






FLA 


0.04 


-0.11 


0.45 


-0.28 


0.64* 








SA 


-0.37 


-0.27 


0.74* 


-0.19 










FCSA 
(<65 |Lim) 


-0.20 


-0.76* 


0.44 












ACSA 
(<65 \im) 


-0.49 


-0.67* 














MeanCT 


0.34 

















^Significance levels: *p<0.01, **p<0.001, ***p<0.0001. 
4. Discussion 



4.1. Minimum cap thickness 

A fibrous cap (FC) with the minimum cap thickness <65 |Lim has been used as the threshold to 
characterize TCFA by histopathology [4], but no causal relationship can be derived for 
prediction of plaque rupture. While largely accepted, the determination of this value may be 
affected by artifacts from histology preparation such as tissue shrinkage, and selection bias 
(this cut-off was derived from plaques that had already ruptured). A recent OCT study on 266 
lesions suggested that 80 \im might be the critical cap thickness that may confer instability to 
a coronary plaque in living patients [17]. However, OCT and pathology studies were based on 
single or multiple individual FC measurements and significant variations can be expected, as 
shown in the present study (Fig. 3). The MAD ofl2.8±8.2|im among the three analysts was 
actually below the resolution of the OCT system. However, the reported incidence of TCFA 
can differ dramatically if 65 jim was chosen as the single threshold (Fig. 3(b)). One particular 
challenge for human analysts is to identify the particular frame containing the cap with the 
minimum thickness. In this study, while the three analysts chose the same frame for most of 
the lesions, the difference in selection of the site to measure the minimum FC thickness could 
be substantial for some lesions (e.g. lesion 2 in Fig. 3(a)). The computer-aided method tested 
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in this study provides semi-automatic segmentation of the whole FC, which eUminates the 
uncertainty associated with manual analysis and image selection and allows more accurate 
determination of the true minimum FC thickness. 




Fig. 7. Two coronary arteries used in the validation study with fibrous caps (FC) rendered in a 
continuous color map indicating the thickness. Both of the two lesions contain TCFA (red 
arrows) and similar minimum cap thickness. However, the one shown in the lower panel 
contains a significantly larger surface area with thin cap as compared to the one in the upper 
panel. 

4.2. Do we know the true FC thickness? 

Xu et al. [25] has found that lipid plaques have higher backscattering coefficient as compared 
to FC. Therefore, the study suggested that the true FC boundary is actually within the signal- 
rich region. If this is true, the FC boundary may precede the steepest decline, and both human 
analysts and our method may overestimate the FC thickness slightly. However, direct search 
for the true boundary as a local peak may not be straightforward for human analysts or for 
algorithms. In clinical images, such peaks are not always present, and even when present, are 
not readily discernible from surrounding speckle noises. In addition, the presence of 
macrophages can also cause very strong peak signals. In comparison, searching the maximum 
gradient and intensity difference to represent the boundary is relatively robust and consistent, 
which is more important clinically than a slight bias. We will perform further validation 
studies to confirm the true FC location. As tissue shrinkage and the difficulty for precise 
registration associated with histology may make it difficult for such validation, a volumetric 
validation method using Cryo-imaging [26] that gives better registration may be used. If 
appropriate, an offset value may be added in the model to compensate for the bias. 
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4.3. Volumetric quantification 



A fundamental limitation of the single number description of the FC thickness is the lack of 
appreciation for the 3-D morphology of FC. While the mechanisms associated with FC 
rupture remain elusive, it is possible that mechanical stability of the cap may not only depend 
on the focal point thickness, but also on the thickness of non-focal regions. Moreover, plaque 
rupture may not necessarily happen at the thinnest part of the FC. These unknown 
morphological and mechanical properties that ultimately lead to plaque rupture can now be 
investigated in future studies with the help of the volumetric quantitative method. 

Among the volumetric metrics, the FC surface area in a certain thickness category reveals 
the component of the FC in a certain thickness range. One may speculate that a FC with a 
larger area of thickness<65 jim may be more vulnerable than a FC with a smaller area of 
thickness<65 jam. The total and categorical FC fractional luminal area is more precise than 
using the number of quadrants to represent the distribution of lipid-rich plaques. The FC 
volume is strongly correlated with the FC surface area. This is not surprising because both 
metrics depends more on the lesion length rather than the lesion thickness due to the large 
frame interval (200 |im) of the current OCT system. We could speculate that vulnerability 
increases with surface area and decreases with thickness, so volume may be less sensitive than 
surface area even though it includes more information. 




^9 




Zoom in 


Auto 




Fig. 8. In this case, the computer algorithm mistakenly identified luminal blood as part of the 
tissue due to the contact between the blood and lumen boundary. In the blood-free region, the 
algorithm accurately segmented the boundary as compared to human analysts. 

4.4. Performance and potential limitations of the method 

The semi-automatic method is as accurate as expert analysts (Fig. 4). It is highly robust 
(global optimum nature of dynamic programming) and fast enough (<ls for single frame) to 
be used in practice. 

One limitation of the semi-automatic method is the inaccurate segmentation of the luminal 
boundary in presence of luminal blood (Fig. 8). In this study, only one lesion had residual 
luminal blood. For validation purposes, no manual correction was performed in this study. 
However, in practice, users could make necessary corrections in such cases. 

The current method is semi-automatic because the circumferential boundaries of a FC 
were defined manually by the analysts. Automated detection of the circumferential boundaries 
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is challenging because validated (or at least reproducible) criteria for the FC circumferential 
boundaries in OCT images have not been defined. 

The metrics provided in this study only characterize the morphology of the FC. However, 
the vulnerability of the FC also depends on inflammation (macrophages) within the FC [27] 
and other factors in addition to FC thickness. For example, Ohayon et al. [28] has 
demonstrated that plaque instability is affected by a combination of cap thickness, necrotic 
core thickness and the arterial remodeling index. The direction of blood flow towards the FC 
may also affect the vulnerability [11]. Microcalcification within the thin fibrous cap has also 
been suggested to cause stress-induced plaque rupture [29]. 

5. Conclusion 

We have demonstrated a computer-aided method that is able to segment the boundaries and 
quantify the 3-D morphology of fibrous caps. The volumetric metrics provide complimentary 
information to the conventional metric (minimum cap thickness). The method is fast, accurate 
and robust, and can be used in future studies for more comprehensive characterization of 
fibrous caps and correlate with clinical presentations to advance our understanding of the 
mechanisms leading to plaque rupture. 
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